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Abstract 

Q_i| The annihilation-creation operators of the harmonic oscillator, the basic and most 

important tools in quantum physics, are generalised to most solvable quantum me- 
chanical systems of single degree of freedom including the so-called 'discrete' quantum 
^3 \ mechanics. They admit exact Heisenberg operator solution. We present unified defi- 

nition of the annihilation-creation operators (o^) as the positive/negative frequency 
parts of the exact Heisenberg operator solution. 

^ . PACS : 03.65.-w, 03.65.Ca, 03.65.Fd, 02.30.1k, 02.30.Gp. 

Introduction 

The annihilation-creation operators are the simplest solution method for quantum mechani- 
cal systems. In this Letter we provide unified definition of annihilation-creation operators for 
most solvable quantum mechanical systems of one degree of freedom. A quantum mechanical 
system is called solved or solvable if the entire energy spectrum {S n } and the corresponding 
eigenvectors {<p n }, 7~L(fin = £ n (fin are known PQ. This is the solution in the Schrodinger picture. 
We will show that they also possess exact Heisenberg operator solutions. The annihilation- 
creation operators are defined as the positive /negative frequency parts of the Heisenberg 
operator solution and they are hermitian conjugate to each other. This method also applies 
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to the 'discrete' quantum mechanical systems, which are deformations of quantum mechanics 
obeying certain difference equations; see full paper [2]. Our results will be translated to those 
of the corresponding orthogonal polynomials. In particular, the exact Heisenberg operator 
solution of the 'sinusoidal coordinate' corresponds to the so-called structure relation [3] for 
the orthogonal polynomials, including the Askey- Wilson and Meixner-Pollaczek polynomials 

m 

Heisenberg Operator Solution 

We will focus on the discrete energy levels, finite or infinite in number, which are non- 
degenerate in one dimension. For the majority of the solvable quantum systems PQ, the n-th 
eigenfunction has the following general structure (f> n (x) = <f)o(x)P n (r](x)), in which 4>o(x) is 
the ground state wavefunction and P n (rj(x)) is an orthogonal polynomial of degree n in a 
real variable rj. 

Our main claim is that this rj(x) undergoes a ' sinusoidal motion' under the given Hamil- 
tonian H, at the classical as well as quantum level. The latter is simply the exact Heisenberg 
operator solution. To be more specific, at the classical level we have 

{H, {H, V }} = -nR (H) - R- X {H). (1) 

The two coefficients Rq and R_i are, in general, polynomials in the Hamiltonian TC. This 
leads to a simple sinusoidal time-evolution: 



» ? (x;t) = 53((-t)Vn!)(ad70 B *7 



n=0 



, sin \tJR (H ) 



Rq{Hq 



+ (rj(x) + R-i(Ho)lBo(Ho)) cos[ty/Ro(Ho)], (2) 

in which a.dli,X = {7i, X} and the subscript means the initial value (at t — 0). The 
corresponding quantum expression is 

[H, [H,Tj] ] = T]R (TC) + [H,Ti] Ri(H) + R-iW, (3) 
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in which R\{Ji) is the quantum effect. The exact Heisenberg operator solution reads 

oo 

e itn r](x)e- itn = ^((rt) n /n!)(ad H) n r] 

n=0 

e ia + (H)t „ict-(H)t 

= [KM*)} JTT, ~ R-l{U)/Ro{H) 

a + {ri) — a^{rt) 

+ („(*) + R-i{H)/MH)) x } a+{n) _ a ^ , W 

in which ad is now a commutator a,d7iX = [H, X], instead of the Poisson bracket. The two 
"frequencies" are 

a± {H) = (RtiH) ± v /R 1 (H) 2 + 4R Q (H))/2, 
a + (H)+a4H) = R\{H), 

a+{H)a-{H) = -R {H). (5) 
If the quantum effects are neglected, i.e. R\ = and 7i — > Ho, we have a + = —a- = 



\/ Ro(Ho), and the above Heisenberg operator solution reduces to the classical one (j2J). When 
the exact operator solution (J3J) is applied to 4> n , the r.h.s. has only three time-dependence, 
e ia±(s n )t an( j a cons t a nt. Thus the l.h.s. can only have two non-vanishing matrix elements 
when sandwiched by <p m , except for the obvious <ft n corresponding to the constant term. In 
accordance with the above general structure of the eigenfunctions, they are <p n ±i', that is 
(4>m\v{ x )\4>n) — 0, for m 7^ n ± 1, n. This imposes the following conditions on the energy 
eigenvalues 

£n+l £"n ^+(^71)5 £fi— 1 &n \&n)- 

These conditions, together with their 'hermitian conjugate' ones, combined with the ground 
state energy £ = determine the entire discrete energy spectrum as shown by Heisenberg 
and Pauli. In this letter we adopt the factorised Hamiltonian: 

U = A j A/2, A4> = 0^ H4> = 0, £ = 0. 

The consistency of the procedure requires that the coefficient of e ia -^* on the r.h.s. should 
vanish when applied to the ground state 0o : 

-[H,r)(x)}<f>o + (v(x)a+{0) - i2_i(O)/a_(O))0 o = 0, 
which is the equation determining the ground state eigenvector O in the Heisenberg picture. 
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Thus we arrive at a dynamical and unified definition of the annihilation-creation op- 
erators: 

e itn r]{x)e- ltn = a^ + \n, n)e ia+{H)t - R_ X (H) / Rq{H) + a^{H, 7])e ia ~ {n)t , (6) 

= (±[H,v(x)} =f ( V (x) + R^{H)/R Q {H))a T {H)) / {a + {H) - a_(W)). (7) 
When acting on the eigenvector <fr n , they read 

a (± Vn(z) = p ^ + (£n - G^lM*) + r~ l{£ i )<Pn{x). (8) 

By using the three-term recursion relation of the orthogonal polynomial P n 

V P n {r]) = A n P n+1 {n) + 5 n P„(r/) + C n P n ^{n) 
on the l.h.s. of ©, we arrive at 

a i+) 4> n = A n <f) n+1 , a (_) 0n = C n <p n _x- (9) 

Based on these relations, it is easy to show that are hermitian conjugate to each other. 
Sometimes it is convenient to introduce a'^ with a different normalisation 

fl /(±) djf Q (±)( a+ ( W ) _ a _( W )) (10) 

= ±[W,77(x)] =F {7i(x)+R-. 1 (H)/R (H))a :f :(H), 

which are no longer hermitian conjugate to each other. 

In the literature there is a quite wide variety of proposed annihilation and creation op- 
erators [Sj. Historically most of these operators are connected to the so-called algebraic 
theory of coherent states, which are defined as the eigenvectors of the annihilation opera- 
tor (AOCS, Annihilation Operator Coherent State). Therefore, for a given potential or a 
quantum Hamiltonian, there could be as many coherent states as the definitions of the an- 
nihilation operators. In our theory, on the contrary, the annihilation-creation operators are 
uniquely determined except for the overall normalisation. In terms of the simple parametri- 
sation ip = Y^=o c n ( t>n{.x) the equation a^tp = Xip, A G C can be solved with the help of 
the formula Q 

00 A n 
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To the best of our knowledge, the 'sinusoidal coordinate' was first introduced in a rather 
broad sense for general (not necessarily solvable) potentials as a useful means for coherent 
state research by Nieto and Simmons jSj. In |2] the necessary and sufficient condition for 
the existence of the 'sinusoidal coordinate' is analysed within the context of ordinary 
quantum mechanics 7i = p 2 /2 + V(x). It turns out that the potential can be expressed in 
terms of the t)(x) and drj/dx 

V(x) = + r { % + c^/{dr]/dx) 2 - rj8. (12) 

The parameters appear in Ro(H) = r^H + 7g , R\{7i) = r± and R-i(H) = r_{TC + r_{ 
and c is the constant of integration. These potentials are all shape invariant j7j. The Kepler 
problems in various coordinates and the Rosen- Morse potential are not contained in ()12j) . 
though they are shape invariant and solvable. 

We give three typical examples. See [2j for more. 

Poschl- Teller potential 

The Hamiltonian of the Poschl- Teller potential, the eigenvalues and the eigenfunctions are 
(0 < x < vr/2): 

Ti = {jp — ig cot x + ihtanx)(p + ig cot x — ihtsmx)/2, 
£ n = 2n{n + g + h), g, h > 0, r](x) = cos 2a;, 

(f) n (x) = (sinx) 9 (cosx)^ Q ' /3) (cos2x), (13) 

in which PlC'^iff) is the Jacobi polynomial and a = g — 1/2, (3 = h — 1/2. The classical 
solution of the initial value problem is (Ti 1 = f Ti + (g + h) 2 /2): 

cos2x(i) = (cos2x(0) + 9 ~, k ) cos[2t v / 2l : F ] 

-p(0)sin2x(0) 1 p_ J - oay , ■ (14) 



2H' 

The corresponding quantum expressions are 

[H, [H, cos 2x\ ] = cos 2a; (8ft' -4) +4{H,cos2x] + 4(a 2 - p 2 ), (15) 



with a±(H) = 2 ± 2^20.'. The annihilation and creation operators are 

J 2 _ Q2 

a'W/2 = a^hVm = ± sin 2x— + cos 2x + - H . (16) 

dx VTR?±1 
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When applied to the eigenvector n as 2£ n + (g + h) 2 = (2n + g + h) 2 , we obtain: 

a /20 "= 2n + a + P ^ (17) 

a'^/Un = 4(n + l)(n + a + /3 + l) 

7 rn 2n + « + /? + 2 rn+1 v ; 

Deformed harmonic oscillator 

The deformed harmonic oscillator is a simplest example of shape invariant 'discrete' quantum 
mechanics. The Hamiltonian of 'discrete' quantum mechanics studied in this letter has the 
following form jHj (with some modification for the Askey- Wilson case): 

H= (^JV{x) e p ^V{xy + ^V{xYe- p ^V{x)-V(x)-V(xy^ I 2. (19) 

The eigenvalue problem for TC, Hep = Sep is a difference equation, instead of a second order 
differential equation. Let us define S±, T± and A by 

S + = e^y^VW, S_ = e- p / 2 ^/7(x), 
T + d ^ f SlS + = ^/V(xje p ^/V(xr, 
T_ = SlS^ = ^/V{xfe- p ^/V{x), 

def 



A = i(S + - S-), A ] = -i{S\ - Si). (20) 

Then the Hamiltonian is factorized 

U = (T+ + T_ - V(x) - V(x)*)/2 

= (S{-Sl)(S + ~S^)/2 = A^A/2. (21) 

The potential function V(x) of the deformed harmonic oscillator is V(x) = a + ix, —oo < 
x < oo, a > 0. As shown in some detail in our previous paper it has an equi-spaced 
spectrum (£ n = n, n = 0, 1, 2, . . .) and the corresponding eigenfunctions are a special case of 
the Meixner-Pollaczek polynomial P„ (x ; | ) jl] , 

0o(x) = \/T(a + ix)T(a — ix), r](x) = x, (22) 

<f> n {x) = Mx)Pn(x), P n (x) = (x ] f ), (23) 

which could be considered as a deformation of the Hermite polynomial. 



6 



The Poisson bracket relations are {TC, x} = —a/ a 2 + x 2 sump, {H, {Tt, x}} = —x, leading 
to the harmonic oscillation, 



x 



(t) = x(0) cos t + \fa 2 + x 2 (0) sinh p(0) sin t, (24) 



which endorses the naming of the deformed harmonic oscillator. The corresponding quantum 
expressions are also simple: [H, x] = —i(T + — T_)/2, [H, [H, x) ] — x, 

e itH xe -WH = x cost + i[n,x] sint = x cost + (T+ - T_)/2sint. (25) 

The annihilation and creation operators are 

a'W = 2a {±) = x±[H,x] = xt *(T+ - T_)/2. (26) 

These operators were also introduced in jH] by a different reasoning from ours. The action 
of the annihilation creation operators on the eigenvectors is 

a'(-)0 n = (n + 2a - l)0 n _ 1; a' (+) 0n = (n + l)0„+i. (27) 

From these it is easy to verify the su(l, 1) commutation relations including the Hamiltonian 
H: 

[H, a'W] = ± a l{± \ [a?(-\ a'«] = 2(W + a). (28) 
The coherent state (jlip . is simply obtained from the formula (J27|) and a'^~' = 2a^: 

^(x) = o (x)f^^pW(x;f), (29) 

which has a concise expression in terms of the hypergeometric function \F± 

^j(x) = <j> Q {x) e 2lX x Fi ( a + 2 * X | -4zA) . (30) 
Askey- Wilson polynomial 

The Askey- Wilson polynomial belongs to the so-called g-scheme of hypergeometric orthog- 
onal polynomials jlj. It has four parameters 01,02,03,04 on top of q (0 < q < 1), and 
is considered as a three-parameter deformation of the Jacobi polynomial. As a dynami- 
cal system, it could be called a deformed Poschl- Teller potential. The quantum-classical 
correspondence has some subtlety because of another 'classical' limit q — > 1. 



The Hamiltonian of the Askey- Wilson polynomial has a bit different form from ()19j) : 

n = (VW)q D VvW+ ^vWq~ D VW)-v(z) - v{zf) / 2, (31) 

with a potential function V(z), z = e lx , < x < n : 

T/ , v ]lj=i(l-a^) n def d .d 

v{z) = (T^kT^) ' D = z Tz = ~ % Tx =p - 

We assume —1 < ai, 02, 03, 04 < 1 and 01020304 < q. The eigenvalues and eigenfunctions are 

m 

£n = (q~ n - 1)(1 - a 1 a 2 a 3 a 4 q n - 1 )/2, 
(f) (x) = 



(z 2 ;g)oo (^ 2 ;g) c 



n i= i(%2 ; g)oo EL^i^ -1 ! <?)oo 
77(2;) = (z + z~ 1 )/2 = cosx, 4> n (x) = (fio(x)P n (cosx), 

P n {v)=Pn{V'^l, a 2,Cl3,a4\q), (32) 

in which p n (r] ; ai, a 2 , 03, a^\q) is the Askey- Wilson polynomial 

The classical sinusoidal motion (J2J) holds with the Hamiltonian (7 = logg): 

H c =y/V c {z)V c {z)*cosh 1 p- (V c (z) + V c (z)*)/2, 
4 

V c (z) = H(l-a jZ )/(l-z 2 f. (33) 

The coefficients in the classical expression are 

R (H c ) = 1 2 (H 2 c + c 1 H c + c 2 ), 
R-i(H c ) = -7 2 (c 3 Hc + c 4 ), 

with ci = 1 + 64, c 2 = (1 — &4) 2 /4, c 3 = (61 + &3)/4, c 4 = (1 — 64) (61 — &3V8. Here we use the 
abbreviation 

4 

, dof \ , dcf \ ^ , def TT 

"1 = a i ' 3 = a J a kai, »4 = ]_]_%■• 

i<i<4 i<j<fc</<4 j=i 

The corresponding quantum expressions are 

Ro{H) = q{q- 1 - l) 2 ((«')* - (1 + g- 1 ) 2 6 4 /4) , 
^(ft) = q ( q ~ l - l) 2 ft', W = H + (l + q-%)/2, 
R-i(H) = -qiq- 1 ~ l) 2 ((6i + q-%)H/A + (1 - q^b^h - 6 3 )/s) . (34) 



The two frequencies are: 

a±(H) = (q- 1 - 1) ((1 - q)H' ± (1 + ?)#Frt) / 2. (35) 
The annihilation-creation operators are: 

a (±) = (±{q- 1 - l)(z~\\ - qz 2 )T + + z(l - g^ 2 )T_)/4 

=F cosx a T (H) ± iL.i(W)a±(W)~ 1 ) / (a+(W) - a-(W)) • (36) 

Their effects on the eigenvectors are: 

( _) _ (1 - q n ) ni< J <fc<4( 1 - "/'*•'/" ' 



2(1 -6 4 g 2 «- 2 )(l -b A q 2n ~ l ) 
1 - b A q n - 1 

2(1 -6 4 g 2n - 1 )(i-M 2riN 



(37) 



The coherent state is 



= o (x)> - = , -P n (cosa;). (39) 

^ (? ; q)n Ui<j<k<^ a j a k ; ?)n 

Conclusions We have shown that most solvable quantum mechanics of one degree of free- 
dom have exact Heisenberg operator solution. The annihilation-creation operators (a^) are 
defined as the positive /negative frequency parts of the exact Heisenberg operator solution. 
These (a^') are hermitian conjugate to each other. This method also applies to the so-called 
'discrete' quantum mechanics whose eigenfunctions are deformations of the classical orthog- 
onal polynomials known as the Askey-scheme of hypergeometric orthogonal polynomials. 
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